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Abstract 

Nonlinear waves have been observed in synchrotrons for 
years but have received little attention in the literature. 
While pathological, these phenomena are worth studying 
on at least two accounts. First, the formation of solitary 
waves may lead to droplet formation that causes significant 
beam halo to develop. It is important to understand the 
conditions under which such behavior may be expected in 
terms of the machine impedance. Secondly, a variety of 
nonlinear processes are likely involved in the normal sat- 
uration of unstable oscillations, leading to the possibility 
that low-level, but potentially broadband fluctuation spec- 
tra may develop. The resulting fluctuation spectra carry 
indirectly the signature of the machine impedance. In this 
work we review a number of observations of nonlinear lon- 
gitudinal waves made in Fermilab accelerators, and make 
a first attempt to develop appropriate theoretical models to 
explain these observations. 

1 INTRODUCTION 

Over the years, nonlinear wave phenomena have received 
scant attention in high energy synchrotrons, in part, be- 
cause of the mathematical difficulty of this subject, but also 
due to the fact that nonlinear wave motion is usually asso- 
ciated with a pathological state of an accelerator that is best 
to be avoided. While this is indeed true for the most violent 
nonlinear effects, a broad class of low-level processes may 
be playing a role in many present machines, and the drive 
for ever higher beam intensities may lead to the widespread 
occurrence of nonlinear wave phenomena. 

In particular, in the case of the dynamical behavior in 
the vicinity of an intense, stored beam, we are interested in 
the formation of beam halo, either as a diffuse cloud, rep- 
resented by a departure from a Gaussian distribution, or as 
droplets which may occur in a type of phase transition at 
beam's edge due to coherent modes. In addition, it is use- 
ful to study the formation of an equilibrium state, if it ex- 
ists, between a broad spectrum of marginally stable modes 
and some weak dissipative mechanisms that can lead to a 
saturated state of low-level turbulence, which, in turn, can 
affect the rate at which the halo population is generated. 
These phenomena can be expected to be most prevalent in 
hadron rings owing to the weakness of the damping mech- 
anisms, and our attention in this paper is focussed on this 
case. 

While these subjects are mathematically complex, a rich 



literature already exists in the field of plasma physics that 
deals with these questions, although the interparticle force 
predominantly considered in this literature is due to space 
charge alone. At the relativistic energies typical in modern 
accelerators, the interaction between particles is dominated 
by wall image currents, i.e. the wakefields, which com- 
plicates the nature of the interaction, but can also lead to 
a wider variety of wave phenomena. It is our aim in this 
work to highlight observations of nonlinear wave phenom- 
ena in high-energy synchrotrons and to point out methods 
of analysis from plasma physics that can be applied to the 
study of these topics. 

The types of wave behavior in beams may be classified 
according to the degree of nonlinearity, in parallel with the 
concepts in plasma physics. In the linear regime, a reso- 
nant mode can be driven resulting in a response at the drive 
frequency which is characteristic of the beam intensity and 
the nature of the wake field, or impedance, of the machine. 
When detected by a suitable pick-up, the driven response 
can shed light on the properties of the wakefields and the 
proximity of the beam to the stability threshold. This so- 
called transfer function method is widely used to study 
accelerator stability. 

If an accelerator is operated just above its stability limit, 
the most unstable mode is driven into exponential growth 
by the wakes, reaching a saturated, though marginally sta- 
ble, state as the beam distribution is altered by the growing 
waves. If the spectrum of unstable modes is sufficiently 
broad, the phase of the perturbation is effectively random, 
and the interaction of waves and particles leads to parti- 
cle diffusion in phase space, known in plasma physics as 
quasi-linear diffusion [Qj. The analog in beams, known as 
the 'overshoot' phenomenon, has been studied al- 
though the applicability of this model is unclear owing to 
the typically narrow unstable spectrum found in many stor- 
age rings. A recent numerical study of this phenomenon 
that shows the complexity of the interaction is found in ref. 

However, particularly in hadron rings where the absence 
of synchrotron damping allows virtually unimpeded mode 
growth, unstable waves can grow to finite amplitudes that 
permit a significant fraction of the beam to become trapped 
in its own wake. The resulting wave motion couples to the 
trapped particles in such a way as to give rise to slowly 
damped oscillations. This phenomenon is known as non- 
linear Landau damping in the plasma physics literature [^j 
(in comparison to linear Landau damping which is part of 



the linear beam response). It is to be expected that where a 
discrete spectrum of unstable waves can occur, as is often 
the case in a synchrotron, that nonlinear Landau damping 
can play an important role. 

At the next higher level of nonlinear interaction, coher- 
ent modes can resonantly interact in a process known as 
the three- wave interaction ftfa. This leads typically to a 
cascade in frequency which, due to the harmonic charac- 
ter of many modes in storage rings, readily occurs and can 
cause a broadening of the original unstable spectrum. This 
phenomenon has been studied in the simple case of longi- 
tudinal oscillations in a coasting beam, [||| and it can be 
expected that similar wave-wave coupling can occur in the 
transverse plane and in bunched beams as well, albeit with 
different resonance conditions. 

If the coherent motion is particularly violent, and suffi- 
ciently dissipative, then the trapped portion of the beam can 
self-extract from the core of the beam distribution, forming 
droplets at the beam edge that can be self-sustaining. These 
are, presumably, a form of solitary wave, or soliton, which 
is perhaps unique to a high-energy synchrotron due to the 
complex character of the wake field. Such solitary waves 
may be a primary producer of halo particles for weakly- 
damped hadron rings. 

In general, we are interested in the final state of these var- 
ious nonlinear interactions: the condition where the coher- 
ent modes reach marginal stability through either a change 
in the beam distribution, or through frequency spreading 
of the spectral distribution. In the latter case, the phonons 
themselves can be thought of as comprising a fluid which 
comes into equilibrium, the details of which depend on the 
inter-phonon interaction. In the plasma physics literature, 
scaling laws for the resulting turbulent fluctuation spectrum 
have been derived ([||] and references contained therein). 
For our purposes, we would like to understand how aspects 
of the machine impedance, and therefore the detailed de- 
sign, contribute to the form of the equilibrium turbulent 
spectrum, if it exists. 

In this paper we review the observations made at Fermi- 
lab [^[ in stored high energy hadron beams and compare 
the observations with numerical simulations. Our exper- 
imental studies, and thus our theoretical work, have been 
focussed on the phenomena in perhaps the simplest of all 
cases, that of longitudinal oscillations of a coasting, or un- 
bundled, beam in a storage ring. As such, the surface of 
this subject has only been scratched, and our aim here is to 
outline the steps that would have to be taken to study any of 
the many other possible situations where nonlinear waves 
can occur. Moreover, we would like to underscore the im- 
portance of understanding turbulence in beams that we feel 
will be playing an increasingly important role as beams are 
more commonly run close to, or even above, their linear 
stability boundaries. 



2.1 Stability in Particle Beams 

In the case of a high-energy stored beam, the growth of co- 
herent wave motion is normally undesireable. Wakefields 
can drive such waves, though the mode growth is counter- 
acted by damping due to the spread in frequencies of the 
individual particles making up the beam, and this damp- 
ing effect was first derived for a plasma, known as Landau 
damping JlCj]. A well known technique for determining 
the linear stability boundary of a beam is to excite driven 
oscillations on the beam and to monitor the amplitude and 
phase of the beam's response, which includes the effects of 
wakefields. This technique, known as a beam transfer func- 
tion, [[!]], yields for longitudinal motion in a coasting beam 
a response of the form 
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where u> s is the harmonic revolution frequency, f is the 
longitudinal particle distribution function, e is the energy 
deviation, k a is the frequency dispersion factor and Z(Q) 
is the machine impedance. This function is directly related 
to the dispersion relation for longitudinal modes given by 



D m (Q) = 1 + Z(Q) 
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Figure 1: Theoretical shift of the beam response due to 
an impedance. The curve centered on the origin (dashed) 
is the response when there is no impedance, and the dis- 
placed curve (dotted) is the response when there is an 



impedance of (Z x , Z y ) = (.05, — .04). 



The magnitude 
.06417, and the 



of the impedance is |M| = <J Z"* 

phase is 9 = —39°. The beam distribution used was Gaus 
sian in energy, and the beam parameters were arbitrary. 



2 REVIEW OF BASIC PHENOMENA 



The stability boundary can be depicted in the impedance 
plane as the curve for which 7m(0) = 0, as shown in Fig.|IJ 



The machine impedance can be extracted from the mea- 
surements as an offset of the centroid of the stability curve, 
provided the beam distribution is known, assumed to be 
Gaussian here. 

2.2 The Three-Wave Interaction 

Weakly nonlinear processes are described using the same 
techniques as in linear stability theory, with the excep- 
tion that a second-order frequency mixing term is in- 
cluded in the description of the dynamics. The effect of 
the frequency-mixing leads to a resonant coupling phe- 
nomenon by which modes at two separate frequencies cou- 
ple to produce a response at a third frequency, a process 
known as three-wave or parametric coupling. The process 
is characterized by selection rules such that 
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corresponding to conservation of energy among the waves. 
A similar condition applies to the mode wavenumbers, cor- 
responding to conservation of momentum. Due to the peri- 
odicity in a ring, this condition can be readily satisfied for 
a large number of normal modes. We have studied the cou- 
pling for longitudinal modes theoretically and found three- 
wave coupling obeys a dispersion relation that couples the 
linear response of harmonic m and m-n through an idler 
mode at harmonic n. 
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where the drive frequency is at fi = noj s and 



Figure 2: Three- wave coupling spectrum for longitudinal 
modes in a coasting beam in the Tevatron, 150 GeV beam. 
Excitation at h = 1000, (47.712 MHz) as shown in the up- 
per graph, led to successive excitation of lower sidebands 
accompanied by low frequency modes, shown in the lower 
graph, which satisfy the selection rule. The lower graph be- 
gins at zero frequency and in each figure the frequency span 
is 2 MHz. The vertical amplitudes are in arbitrary units but 
the scales are logrithmic. 
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and V is the drive amplitude, I is the beam current, r/ is 
the slip factor, #- is the fractional energy spread and E a is 
the beam energy. 

The implication of Eq. Q is that three-wave coupling is 
most likely near the stability threshold for any of the modes 
involved. The selection rule Eq. [3] leads to a single-sided 
coupling, which was observed experimentally, as shown in 
Fig- I 

An interesting issue to investigate is how the power in 
the excited modes varies in time, especially in the presence 



of damping. Experimental observations indicate that a very 
regular cascade toward lower frequencies takes place, evi- 
dently due to successive three-wave coupling events. This 
behavior is typical for a dissipative system with sufficiently 
high mode density, and may be described by the following 
system of equations for the mode amplitudes. 

— 7t, — = s m e i ^V mri ) : A n Ak (5) 
at 

where the matrix element of the interaction has been sym- 
bolized as V mn k, and is defined by as the following, 
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Figure 3: Power versus time at h = 105 in response to a 
0.5 msec drive pulse. The placement and duration of the 
drive pulse has been drawn in for reference. An impulse 
excitation leads to slowly decaying amplitude oscillations 
as trapped particles exchange energy with the wakefield. 
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It is worthwhile to note that as the multiplicity of modes 
becomes sufficiently dense, the coupling between waves 
governed by Eq. BJcan lead to a solitary wave phenomenon 
QlTp, and this subject will be described further in a later 
section. In the above mentioned work, only the interac- 
tion of longitudinal modes has been considered. It is also 
reasonable to expect that transverse modes can be coupled, 
especially where nonlinearities can play an important role, 
such as in the beam-beam interaction. This should be a 
fruitful area for further study. 

2.3 Nonlinear Landau Damping 

Sufficiently large wakefields disturb an initially smooth 
particle distribution by trapping particles within the po- 
tential wells of the waves generated. The particle motion 
decoheres with a time constant that is significantly longer 
than the inverse frequency spread, or linear Landau damp- 
ing time. The trapped particles undergo synchrotron os- 
cillations in the self-generated potential wells, alternately 
exchanging incremental energy with the wakefields. The 
combination of energy dispersion of the particles and the 
nonlinearity of the voltage waveform eventually causes 
phase mixing of the coherent motion. This nonlinear damp- 
ing process is called nonlinear Landau damping, and was 
first studied in plasma physics. [Bp, JT^ ] - JT^|. 

Experiments were carried out in the Fermilab Main Ring 
which clearly showed the signature of nonlinear Landau 



damping. In these studies, a short pulse of rf power was ap- 
plied to the beam using an rf cavity at h=106 (5.03 MHz). 
The resulting response showed a characteristic response 
whose envelope decayed not exponentially but in an oscil- 
latory manner, as shown in |[ This behavior is attributed 
to the exchange of energy between trapped particles and 
waves as described above. Both analytic Jig], [|l2|l, [15|, 
and numerical work [Q, on nonlinear Landau damping 
has been carried out which descibes the behavior we have 
observed. This result will also be discussed further in a 
later section on simulations. 

It has also been pointed out [ [To] ] that the advent of parti- 
cle bunching is accompanied by the appearance of coherent 
power in higher harmonics of the fundamental frequency 
of the wakefield as the trapped particle bunches compress 
within the potential wells. Such behavior has indeed been 
observed in the experiments described above Q. This 
compression of the bunch length is essentially wavefront 
steepening which is a prerequisite for the formation of soli- 
tons. 

2.4 Solitary Waves. 

The formation of solitons in a beam is of interest, since soli- 
tons may well be the vehicle which carries coherent energy 
in a highly turbulent state that might occur in a beam with 
weak damping. A vast literature on solitons in various me- 
dia exists [ |l7| | - pojl, though little effort has been given to 
this subject in ultra-relativistic beams. In particular, soli- 
tons in plasmas have been studied extensively, [[[80, [20|, 



which depend on the particular nonlinearity introduced by 
the Coulomb force, i.e. space charge. A similar space 
charge limit was studied for a coasting beam [|TJ, p2|], and 
for a resonator impedance, [|23|], leading to the possibility 
that solitons may exist in high energy beams under certain 
conditions. 

Since a wakefield force is fundamentally more complex 
than the space charge force, it can be assumed that the char- 
acteristics of a soliton will also be unique to the case of a 
high-energy beam. In particular, it is interesting to know 
what the impact of the wakefield dissipation has on the soli- 
tary wave behavior. Results from a variety of beam experi- 
ments suggest that long-lived solitary structures may form 
and extract themselves from the core of the beam over long 
times [|4]], [j25|]. In other work, transient solitary waves 
seem to appear [^J. In all cases, solitons may be viewed 
as phase-space droplets that appear in the beam, and under 
some conditions, give rise to a phase-transition and clumpy 
halo formation. From this point of view, it is valuable to 
understand their dynamics. 

To this end we sketch here the results of an analytic study 
of longitudinal solitons on a coasting beam due to a general 
resonator impedance. This reperesents the simplest possi- 
ble scenario for understanding such phenomena, will illus- 
trate the mathematical procedure and serves as the starting 
point for more complex situations. For the reader's sake 
we note that many steps have been omitted from the fol- 
lowing derivation for reasons of space. Full details will be 



given in forthcoming work |]26[]. The model equations for 
the dynamics are given by the following system 
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where 



I(6;T)= I dvvf(8,v;T), 
where / (6, v; T) is the longitudinal distribution function, 



V is the voltage on a resonator of 7 
being the resonator frequency, and I is the instantaneous 
beam current. Time t has been normalized as T = ui s t. 
Furthermore v — — 4x 
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angular velocity of a beam particle and 
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where R is the resonator shunt impedance. Using standard 
moment techniques [Q] on the above equations, we may 
pass over to the hydrodynamic picture of longitudinal beam 
motion and start from the system of gas-dynamic equations 
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where p and u are the density and the mean velocity mo- 
ments of the distribution, respectively, (the variables p and 
V have been appropriately scaled) and A = p M- (Po = 
constant is the equilibrium beam density.) Using a renor- 
malization group approach [p^], [|2~7j], we may derive a set 
of amplitude equations for the rescaled beam density R a , 
the current velocity u and the mode envelope function E. 

Before proceeding, we would like to examine the stabil- 
ity problem of stationary waves in this system, which can 
be done without a formal solution of the amplitude equa- 
tions. The approach, introduced by Sagdeev JT^], is to look 
for forms of the nonlinear equations which correspond to 
harmonic motion in an effective potential well. Such states, 
if they exist, are conjectured to be allowed solitary waves 
in the nonlinear system. To this end let us write down the 
full system of amplitude equations, which after appropriate 
scaling reads as 
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are the new (rescaled) dependent variables. The coeffi- 
cients entering the above expressions are specified as fol- 
lows: 
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where a v is the normalized beam energy spread. The new 
independent variables (time r and azimuthal position 9) 
are given by 
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Moreover in the above set of equations the following nota- 
tions have been adopted 



= or — A 



222 

= - 7 , 



In order to proceed, we have to further assume that the res- 
onator is weakly damped, namely the high-Q case. For this 
case (7 = 0) the ansatz 
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leads to a system of differential equations for p, v and ip 
admitting the following integrals of motion 
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Figure 4: Pseudo-potential based on the nonlinear wave 
equations for a coasting beam. A minimum in this poten- 
tial indicates the possibility for solitary waves (cavitons) to 
form. 
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These integrals of motion suggest that the stability of sta- 
tionary waves can be equivalently described in terms of 
motion of a single particle in a (pseudo)- potential well. 
Indeed, the function 
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provided p is expressed in terms of A from the second in- 
tegral, comprises a pseudo-potential function. In Fig. ^| we 
show U (A) for the simplest case of constant current ve- 
locity v — v (Ci — 0) . We note that a minimum in this 
pseudo-potential corresponds to solitary waves that can be 
effectively "trapped" in this potential well. 

In the following, we can proceed to find approximate 
closed-form solutions of these nonlinear equations, which 
will allow us to explicitly find the time behavior of the 
solitary waves in the presence of dissipation. Eliminating 
of the current velocity from the complete set of amplitude 
equations and expressing p in terms of \-ip\ 2 
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we finally arrive at the damped nonlinear Schrodinger 
equation 

■ d ^ , ■ hi 
i— — h rybw = 

OT 



2ij\ d 2 tp a-ydij) 2 
1 7TT H 7^ W V>, 

U> n I OX* UJ n OX 




(13) 



Figure 5: Voltage amplitude evolution of the solitary wave 
due to a resonator impedance. In the frame of the wave, 
the amplitude persists for long times but eventually damps 
as the soliton decelerates away from the beam core, and 
hence, the resonator's resonant frequency. 

where 

Eq. [□] admits closed form solutions that indicate solitary 
waves can exist, but due to the dissipation in the model, 
eventually disappear after initial generation. We interpret 
this behavior as a gradual shrinking of the potential well 
that occurs when the trapped particles have decelerated suf- 
ficiently from the the resonantor frequency. The results 
for the voltage amplitude and soliton (caviton) density are 
shown in Figs. |] and ^ respectively. We note that a similar 
equation has been derived [|2[ from an entirely different 
perspective. It is also worth noting that Eq. [j~3| is a special 
case of the complex, cubic Ginzburg-Landau equation 1 29 1, 
widely used to study various pattern formation phenomena 
and coherent structures. 

2.5 Turbulence. 

The study of turbulence in beams is valuable primarily be- 
cause it may be a universal phenomenon, at least at low 
levels, which plays a role in determining the limiting phase- 
space density in any machine. The effect has likely been 
small in machines well below their stability thresholds, 
however, as intensities have been pushed closer to stability 
limits, nonlinear wave interactions can occur which lead to 
a marginally stable equilibrium. A first attempt at deter- 
mining the fluctuation spectrum for a beam was obtained 
by considering the equilibrium state of a Gaussian beam 
[p0|. The resulting spectrum was related to the linear di- 
electric function and showed that the fluctuation density 
would be strongly peaked for cases near the linear stabil- 
ity limit. It is our conjecture that such a situation may have 




Figure 6: Caviton density evolution. A density depression 
associated with the solitary wave decreases in amlitude at 
long times as the wave amplitude itself decreases. 

occurred in the Fermilab Tevatron during recent attempts to 
realize stochastic cooling of bunched beams pip. A broad, 
stationary spectrum of fluctuations was observed at many 
times the expected Schottky, or shot noise, levels. The 
harmonic generation observed in the Fermilab Main Ring 
mentioned above is consistent with these observations as 
well. It is our aim in this paper to outline a theoretical ap- 
proach to understand the formation of an equilibrium spec- 
trum and to understand the spectral amplitude dependence 
on the character of the machine impedance. 

The approach taken is to develop a statistical description 
of fluctuations for an ensemble of coupled modes. This 
may be viewed as a development of the amplitude equa- 
tions associated with coupled modes, as in Eq. | The 
interaction between modes may be three-wave, as given in 
H or higher-order. In this work, we outline a general proce- 
dure, but keep only interactions up to the three-wave level. 
The result will be a scaling law for the envelope of the fluc- 
tuation spectrum. 

The starting point for our analysis is the system of equa- 
tions for the fluctuation 6 N of the microscopic phase space 
density and the fluctuation 8V of the voltage ||] 
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written in the variables a = 6 — uj s t and v = k e/ui s . 
Fourier transforming the above equations and using the 



concept of slowly varying amplitude of weakly nonlinear 
waves one obtains the following equation 
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are the group velocity of waves and the damping factor re- 
spectively. For the slowly evolving part A k of the wave 
amplitude SV k it is straightforward to derive the equation 
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Averaging of the equation for A k yields the kinetic equa- 
tion for waves 
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Dimensional analysis of the kinetic equation for waves 
gives the following fluctuation spectrum law of Kol- 
mogorov type 



const * k 7 / 3 . 
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We note that this power law spectrum, which in this case 
is due to a resonator, would indeed lead to the type of broad 
fluctuation spectrum seen in experiments. However, at this 
time, a detailed study of the scaling of the observed spec- 
trum has not been made. An extension of the above work 
to consider bunched beams, and other types of machine 
impedance would be very worthwhile. 

3 SIMULATIONS 

In this work, we are interested in demonstrating examples 
of the nonlinear wave behavior we have described analyt- 
ically. No attempt is made to closely model a real device, 
though this could be done with the building blocks we pro- 
vide here. We shall concentrate on the longitudinal plane, 
as explained above, and carry out simulations exculsively 
on a coasting, (unbunched) beam, for simplicity. We adopt 
the resonator model described above and follow approx- 
imately the procedure adopted in early simulation work 
[53]. The particle evolution equations are given by 



de 



= 2ireuj s V(t) 



de 



= -m 



dV UJ r R Ul r 



(18) 



dl\ uj r C 



dt 



Figure 7: Simulation after 500 time steps for the model 
problem listed in Table 1. The upper portion represents 
e — 9 phase-space. The middle curve is the cavity volt- 
age as applied to different portions of the beam and the 
bottom curve is the projection of phase-space on the theta- 
axis, representing the instantaneous current. The resonator 
wakefields have caused bunching which shows significant 
particle trapping and consequent wave overturning. The 
particle current shows strong local intensification of the 
beam density corresponding to the solitons. 



where e is the energy deviation from the synchronous en- 
ergy and V is the voltage induced in a resonator with shunt 
impedance R, resonant frequency u r and quality factor Q. 
I is the instantaneous current given as the projection of 
phase-space onto the 9 axis 



1=^1 /* 
2vr ' 



The energy dispersion 77 is assumed to be a constant and 
the particle distribution is advanced each time step accord- 
ing to Eq. Then the current is computed, followed by 
the voltage on the resonator using Eq. Il8[ Results for the 
case of a coasting beam are shown in Figure - |[ The 
simulation parameters are given in Table [l|. 

In Fig. [7], the phase-space distribution is shown, initially 
assumed to be uniform in 9 and Gaussian in energy. Af- 
ter 500 time steps, the resonator has developed a sinusoidal 
voltage from the initial noise level (due to the finite num- 
ber of particles) which has succeeded in bunching a large 
fraction of the beam and synchrotron motion in the result- 
ing potential well is taking place. The synchrotron motion 
can also be thought of as synonymous to the wave-breaking 



R 



-V 



process which has been described in plasma physics [12|. 
As time proceeds, Fig. || shows a deceleration of the 
trapped particles from the core of the beam, and the decel- 
erated portion remains well-organized, and even intensifies 
as its length is foreshortened. The voltage in the resonator 
then becomes phase-locked to the 'droplets' and the voltage 
amplitude oscillates as they move in and out of phase with 






Figure 8: Phase-space after 1000 time steps. Portions of 
the originally trapped particle population have decelerated 
from the core due to the finite resistivity of the wakefields. 
The solitons remain, however, well-organized as they expe- 
rience deceleration. Both the cavity voltage and the density 
perturbations have decreased, but persist for long times. 
The diagrams have the same meaning as in the previous 
figure. 

the remaining coherent structure in the beam's core. We 
note that the droplets thus formed bear the characteristics 
of the solitons discussed in the previous section, remaining 
self-organized for long times. 

As the dissipation in the impedance continues to decel- 
erate the solitons Fig. ^, the resonator voltage drops due to 
the high Q value, or narrow bandwidth, assumed. This, in 
turn, reduces the deceleration rate and the depth of the po- 
tential wells that can sustain the solitons. As such, a steady 
state can be reached where the remaining trapped particles 
reach a stable equilibrium outside the beam, in accord with 
the analytic model for solitary waves. The envelope of the 
cavity voltage and the mean energy deviation of the solitons 
are shown in Fig. [l(| indicating deceleration of the trapped 
particles. A steady state is eventually reached, though not 
shown, where the solitons have moved sufficiently off res- 
onance that the deceleration ceases. 

We show the final results for a low-Q cavity in Fig. [TT| . 
These are qualitatively the same as in the previous case, 
but the structure of the solitons has taken on a decidedly 
random character. This is evidently due to the fact that 
the fractional contribution of noise to the cavity voltage is 
larger, owing to the wider bandwidth of the cavity, result- 
ing in a more random distribution of potential well sizes. 
Droplet, or soliton formation, still occurs, but the resulting 
fluctuation spectrum is broader. The onset of the solitary 
waves can be viewed as a phase transition at the beam's 
edge produced by the resonator's wakefield. This is the 
case, we believe, that is most frequently encountered in ac- 
tual machines. 

We note that the fluctuation spectrum associated with the 



Figure 9: Phase-space after 2000 time steps. The beam 
core is now largely decohered, however, the solitons con- 
tinue to decelerate slowly and maintain the remaining cav- 
ity voltage well off resonance. The diagrams have the same 
meaning as in the previous figure. 



Table 1 : Model Simulation Parameters 



Parameter 


Value 


Resonator Impedance 


100 Ohms 


Slip Factor 77 


.001 


Resonator Quality Factor 


10. 


Typical Beam Energy Spread a £ /Eo 


l-4x 10- 3 


Number of Particles 


10000 



above distribution is due to the ensemble of strongly non- 
linear waves and is likely beyond the realm of the three- 
wave interaction described in the scaling law in the previ- 
ous section. An interesting study to be carried out is the 
experimental and theoretical determination of the spectral 
shape in a machine whose impedance is well-known. 

4 CONCLUSIONS 

In this work we have attempted to outline various levels of 
nonlinearity in coherent interactions in high-energy beams. 
Besides the general academic interest of nonlinear dynam- 
ics, for which high-energy beams provide an excellent test- 
ing ground, there are at least two areas where the study of 
nonlinear waves can find application in accelerator physics. 

The first is the study of halo formation in which the non- 
linear evolution of coherent fluctuations can lead to soli- 
ton, or droplet, formation, as described in previous sec- 
tions. While there is suggestive experimental evidence that 
such states can occur, there has been little detailed study 
of this phenomenon, and we assert that there is much to 
be learned about machine wakefields through the study 
of solitary waves and their interactions. Specifically, we 
have only considered the simplest case, that of longitudi- 
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Figure 10: Cavity voltage envelope and mean energy devi- 
ation. The upper curve is the envelope of the cavity volt- 
age, showing the exchange of energy between the soliton 
and the wakefield. The lower curve is the mean energy of 
the perturbation, which descends away from the beam, but 
eventually reaches an equilibrium (not shown) where the 
deceleration ceases. 

nal waves in an unbunched beam, and there are many other 
cases of interest in high-energy accelerators. 

The second application is the study of non-equilibrium 
fluctuations driven by wake fields. Nonlinear mode-mode 
coupling permits a frequency cascade, both toward lower 
and higher frequencies via separate processes. The pho- 
ton distribution that results is an equilibrium between the 
nonlinear interactions producing the cascade and weak dis- 
sipative mechanisms. These mechanisms are assumed to 
be related to the broadband impedance of the machine, 
though other mechanisms may also be responsible. We 
have carried out a model calculation for a specific form of 
impedance that yields a specific scaling law for the turbu- 
lent spectrum. A number of assumptions come into play in 
the development of this model and the situation is ripe for 
careful experimental testing. The benefit of this study is the 
understanding of the significance of low-level turbulence in 
the limiting parameters of a given accelerator. 
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